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some understanding on the way the fluctuation-dissipation theorem is violated out-of- 
equilibrium. 
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1. Introduction 

The knowledge about out-of-equilibrium processes is far from being as advanced as 
for systems at equilibrium. In particular, the fluctuation-dissipation theorem (FDT) 
which holds at equilibrium is known to be violated out-of-equilibrium. This theorem 
states that at equilibrium the response function Req{t, s) at time t to an infinitesimal 
field branched to the system at time s is related to the time-derivative of the two-time 
autocorrelation function Ceq(t, s): 

Req{t,s) = /3^^^^. (1) 

In the Ising case, the response function reads R{t, s) = "^j^'^*^^^ and the correlation 
function C(t, s) = (crj(t)crj(s)). Based on a mean-field study of spin-glasses, Cughandolo 
et al 1^ have conjectured that for asymptotically large times the FDT has to 
be corrected by a multiplicative factor X{t, s) which moreover depends on time only 
through the correlation function: 

R(t,s)=f3XiCit,s))^^^. (2) 

Exact results have been obtained for the ferromagnetic Ising chain ^ that confirm this 
conjecture. The violation ratio X{t,s) has been also numerically computed for many 
systems: 2d and 3d- Ising ferromagnets 3d Edwards- Anderson model 0, 3d and 4d 
Gaussian Ising spin glasses [0, 2d Ising ferromagnet with dipolar interactions [§[,... 
The response function has been exactly calculated in the case of the Ising chain 
but the result is not written in terms of thermodynamic averages making difficult any 
conjecture on the way the FDT is violated. 

We present an analytic study of the Glauber dynamics of the Ising model with 
arbitrary exchange couplings. The response function to an infinitesimal magnetic field 
is first calculated at equilibrium and the compatibility with the fiuctuation-dissipation 
theorem is verified. In the last section, the response function is derived out-of- 
equilibrium and the expression is compared to the equilibrium equivalent. 

2. Model, dynamics and general response function 
2.1. The model and its dynamics 

We consider a classical Ising model whose degrees of freedom are scalar variables ctj = ±1 
located on the nodes of a dimensional lattice. Let us denote ^({0"}, t) the probability 
to observe the system in the state {a} at time t. The dynamics of the system is supposed 
to be governed by the master equation: 



Y: [piW'}, t)W{{a'} {a}) - p{{a}, t)Wi{a} ^ {a'})] (3) 
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where iy({cr} — > {cr'}) is the transition rate from the state {a} to the state {a'} per 
unit of time. The normation is written I]{o-'} ^({^} ~^ W}) = 1- The transition rates 
are defined by the stationarity condition 

E [Pe,{W'})W{{a'} {a}) - p.,{{a})W{{a} ^ {a'})] = (4) 
W'} 

where peq{{(^}) is the equihbrium probabihty distribution which reads for the Ising 
model with general exchange- couplings: 

Pc,({a}) = le-/^«(M) = 1 -^S^.<^^-'^'='^'. (5) 
Equation (^) is satisfied when the detailed balance holds: 

Peq({a'})W^({a'} ^ {a}) = p.,{{a})W {{a} - {a'}). (6) 

This last unnecessary but sufficient condition is fulfilled by the heat-bath single-spin flip 
dynamics defined by the following transition rates: 



Wi{a} {a'}) 



l^k 



J2a Peq{{(yi}l^k,Crk = ) 



:±1 



(7) 



Only the single-spin flip ak cr'k is allowed. The product of Kronecker deltas ensures 
that all other spins are not modified during the transition. The spin ak takes after the 
transition the new value o"^ chosen according to the equilibrium probability distribution 
pcq{{a}). All spins are successively updated during the evolution of the system. These 
transition rates reduce to Glauber's ones M in the case of the Ising chain. 



2.2. The response function 

A magnetic field hi is coupled to the spin cTj between the times s and s + As where 
As is supposed to be small. Between these instants, the transition rates are changed 
to W/i({o"} — ^ W}) satisfying equation with an additional term (3hiai in the 
Hamiltonian of the equilibrium probability distribution (O): 



l^k 



(8) 



Y,(T=±1 6 

The transition rates are identical to the case /ij = apart from the single-spin flip 
involving the spin cTj. Using the Bayes relation and the master equation @, the average 
of the spin ai at time t can be expanded to lowest order in As under the following form: 

= E ^ip(W,t|M,s + As)p({(7'},s + As) 
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= E (rip{{a},t\{a'},s + As)[{l-As)p{{a'},s) 
WhW'} 

+ Asj:pi{a"},s)W,i{a"}^{a'})]. (9) 
W"} 

Wh being the only quantity depending on the magnetic field in equation (|), the first 
term disappears after a derivative with respect to the magnetic field. In the last term, 
the single-spin fiip involving the spin cTj, i.e. the one coupled to the magnetic field, is 
the only one that gives a non-vanishing contribution. The integrated response function 



between the time s and s + As is 
function: 



dhi 



and the limit As leads to the response 



1 



lim , 

As^O As 



d{a,{t)) 



dhi 



2 j-s+As 

= lim — / Rii{t,u)du = Rii{t,s). (10) 

According to the definition of the response function ([l0|), the calculation can be limited 
to the lowest order in As. As a consequence, the conditional probability in equation (P) 
can be replaced by p{{(j}, s). 



3. Equilibrium fluctuation-dissipation relation 

3.1. The response function 

If the system is in thermal equilibrium at time s, the probability distribution p({(t}, s) 
has to be replaced in equation (^) by peq{{cr})- Then the sum over {a"} can be 
calculated: 



E Pe,{W"})W,{{a"} {a'}) 

{a"} 



nSa",a' 
1 ' 7 



E. , . Jaa'a'. —hia'. 



Pcq({^T'})e 



cosh Jijcr'j 



E . , . Jiwo-'^'cr' —hi 



J3) 



(11) 



cosh (^(3 Y^j^i Jij<^j - hi 

The derivative with respect to hi followed by the limit hi ^ gives (3peq{{a'})[ai — 
tanh {P Jij'^j)]- When replaced in equation (||), one obtains the following integrated 
response function 

^d{a,{t)y 



dhi 



(3As 



hi^O 



{ai{t)ai{s)) - {ai{t) tanh (^p J2 Jij^ji^))) 



and according to equation (|TOD the response function reads 

Rii{t,s) = /3 {a,{t)a,{s)) - (ai(t) tanh (/?^ Ji,a,(s))) 



(12) 



(13) 
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3.2. Equivalence with the usual expression 

The equilibrium fluctuation-dissipation theorem states that the response function is 
related to the time-derivative of the correlation function. We will now show that 
the expression ([13D is compatible with this statement. At equilibrium, the correlation 
function is invariant under a time-translation. As a consequence, one has the relation 
dCii{t,s) _ dCii{t,s) 

dt ~ ds ' ^ ' 

Moreover, it can be shown that the conditional probability p{{a},t\{a'}, s) satisfies the 
master equation for the time-variable t and a similar equation for s: 

(l - ^) p{W},t\{a'},s) = EpiW}^t\W"}^^mW} - K})- (15) 

The time- derivative of the correlation function is then 
dCii{t,s) ^ ,dp{{a},t\{a'},s) 



dt 



dt 



-pcqiW'}) 



(16) 



> dp{{a},t\{a'},s) , 

- ^i^i Peq{{(^ \) 



= -Cu{t,s) + E cr.a'MW},t\W"},s)W{{a'} ^ {a"})p,^{{a'}). 

W,{<7'},{<7"} 

The sum over {a"} can be performed on the lines of equation (p!T| ) for the Ising-Glauber 
model. One obtains 



= peq({a"}) tanh E ■J^J^J j ■ (17) 

Replacing in (|1^) and identifying with ([I3D, one gets the expression of the equilibrium 
fluctuation-dissipation theorem for the Ising-Glauber model: 

Rdt, s) = = (3 {a,{t)a,{s)) - (a,(t) tanh {f3Y.J^M^))) 



4. General fluctuation-dissipation relation 

Calculations have been made up to now within the hypothesis that the system is at 
equilibrium at time s. We will now derive a relation valid even if the system is far from 
equilibrium. We start with the expansion @ of the average {ai{t)). Wh being the only 
quantity depending on the magnetic field, the integrated response function between the 
time s and s + As is to lowest order in As 



^{a^{t)) 
dhi 



AsY. a,p(M,t|{a'},s) 



hi^O {<T},{cr'}, 
{a"} 



dhi 



p{{a"},s). (19) 
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The derivative of the transition rate Wh defined by equation (§) is easily taken and 
reads 



dWh 
dhi 



PW{{a"} {a'}) 



a[ — tanh (3 ^ Jij(y'j 



(20) 



Combining (|TU|) and (^^), the integrated response function is rewritten as 
^^""'^^^^ E a.a',p{{a},t\{a'},s)W{{a"}^{a'})p{{a"},s) 



dhi 



Asp E cT,p{W},t\{a'},s)W{{a"}^{a'}) 



(21) 



{""} 

X tanh ^ Jija'^ piW"}, s). 

In the first term of (|2T|), the sum over {cr"} can be performed using the master equation 
(0) and the expression reorganised using the fact that after the hmit hi 0, the 
conditional probability becomes invariant under a time-translation and depends only on 
the difference t — s. One obtains 

Y: ay,p{{a},t\W},s)W{{a"} - {a'})p{{a"}, s) 



E a,alp{{a},t\W},s) [l + ^] p{{cr'},s) 

o"},|o"'} ^ ' 



W},W'} 

Cu{t,s) + ^^^- E ^^^'^^iW},t\W'},s)p{{a'},s) 

W},W'} 



n (+ \ ^ 9Cii{t,s) dCii{t,s) 
'-'iiy'C, S) -\ h 



3s ' at ■ (^2' 

Similarly, in the second term of equation (|2T|) the sum over {a"} can be performed using 
the master equation (^ and leads to 

Y: a,p(M,t|{a'},.)tanh f /^E^^.^; ) Wi{a"} K})p(K},.) 



E a.p(M,t|{a'},.)tanh(/3E^..^;j + 

) E ^.p(W,t|K},s)tanh(/5E'^^.«^^| (23) 



ds dt 



Putting together equations (|2^ ) and ( p3D into (|21|), one can then evaluate the response 

(24) 



function ([T0|): 



Cji(t, s) - ((Tj(t) tanh ( P E Jij^^ji^) ) ) 
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The hyperbolic tangent is the equihbrium value of a single spin cTj in the Weiss field 
created by all other spins at time s. Moreover, the hyperbolic tangent being an odd 
function of its argument and since af = 1, one can rewrite (E^ in the following form: 




Cii{t, s) - {ai{t)ai{s) tanh 
Hj^i Jij^i{.s)(^j{s) is the energy at site i. 



(25) 



5. Conclusion 

The expression (|2^ (or (pSf )) of the response function makes apparent the way the 
FDT is violated. At equilibrium, all averages are invariant under a time-translation, 
i.e. depend only on t — s so the derivatives over s and t cancel and at equilibrium the 
average in equation (^) is equal to the time-derivative of the correlation function as 
shown in section 3.2. The expression (p^) is very general. Since the response function 
is expressed as thermodynamic averages, it does not depend on the initial conditions. 
Moreover, it does not depend on the dimension of space or whether or not, the lattice 
is regular. It applies to any configuration of exchange couplings so for a disordered 
system, the average response function is just the average of expression (|2^ ) over all 
couplings configurations. Equations (p2D and ( p3D suggest that the presence of the 



term 1 + ^ + ^ in the expression of the response function may be a general feature 
of the response function out-of-equilibrium. Exact results for other models are highly 
desirable. The last term of (|2^) whose equilibrium limit is the time-derivative of the 
correlation function may also be evaluated perturbatively for certain models. 

This work may also help to improve numerical simulations. Up to now, the FDT has 
been numerically checked by calculating the time-integrated response function. One can 
note that equation ( |T9|) can be used to calculate the response function directly during 



a Monte Carlo simulation. 
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